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f""') ■ Abstract 

O^l A complete model of baryon production in an expanding, primor- 
■ dially symmetric hot Universe is constructed in the framework of general- 
ly ' relativistic kinetic theory. In this model specific model for a baryon is 
. calculated and graphs of the value dependence are constructed. 

t> 

1 Transformation to the dimentionless variables 

O In the previous papers of the authors [H , [2] in terms of the developed relativistic 

kinetic theory of baryons production^ an expression for final concentration of 
?-h , baryons in hot Universe was obtained (23) [2J, using which along with the dense 

, entropy equation (26) [2], we get the sought- for expression for final dense of 

specific entropy accounted for one baryon 

S S = ^P = / exp | - / WW I G(t) d t, (1) 

' where 

*(t) = ^| / [f -!>,/;, ;,„/?. i2i 

O ' 



\ t 





oo 



G(t) = 4j / P 2 $<5/dP; (3) 

X. 

derivatives by time are denoted by a dot t, Sf - is a deviation of the boson 
distribution function from equilibrium in a symmetric plasma A = 

t 

Sf(F, t) = - e -*( p >*) / e*( p >*')/ (0; P, t')dt', (4) 







and the following designation is used 



. . 1 f a(t')P (F,t')dt' 
to J ^/a 2 (t')+F 2 /m 2 x 



1 As well as the other similar particles arising as a result of thermodynamic equilibrium 
violation and spontaneous CP invariance violation. 
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According to the formulae (28)-(35)[2] we come from the temporal variable t 
and the impulse one P to the dimcntionless temporal variable 77 and the impulse 
variable £ 

t = r ?y. (6) 

In doing so we choose the normalization of scale factor for a completely ultra- 
relativistic stage of universe expansion, such one that 

a(t) = y/wj, (7) 
where To is the supermassive boson decay time in its own frame of reference 

tq = — ^(m X a) . (8) 

and introduce a dimentionless parameter a 

m x _ rnxy/ro _ X\/ m x 



T(t ) To JctTo 



(9) 



The dimentionless impulse variable £ let us introduce with the help of the rela- 
tion (Ref. [2]): 

P = m x sJToi => i = — ^— (10) 
mxy/To 

so that 

E 



- = yJm* x +Fya*(t)T = a^Te (11) 
and the equilibrium function of the supermassive bosons is equal to (Ref. (13) [2 ) 

e °V+? — 1 

As the done investigations have shown the final results of the investigations are 
very weakly sensitive to statistic factors consideration by the factors $(?y, £) cal- 
culating, whereas they are very sensitive to the static factors consideration at the 
other stages of calculation. Therefore further on we shall calculate the function 
77, £ in the Boltzmann approximation , while at the other stages of calculation we 
shall hold the statistic factors. Then in the Boltzmann approximation f3 ~ 1/2, 
and we obtain for the function $: 

and - 



m v) = \ [vvV^Te - e m ^ + f TF j . (i 4) 

In consequence of (II 3[) G(?7, £) is a nonnegative monotonically increasing func- 
tion 

* MS0; ^Ma-^,., ( 15 , 
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with 

lim $(77, £) = 0. (16) 

77— >o 

Note that the rather large expression for the function <!>(£, 77), obtained in [2J, 
can be reduced to the given above after simple transformations. 

Thus coming to the new variables we get an expression for specific entropy 

00 00 
v 

where 

00 

9(7?) = J V(t)dt. 
t 

As far as $ > 0, then ^(t) > 0, thus 0(77) is a nonnegative monotonically 
decreasing function 

In this case 

n 

Sf(v,0 = e-*<">« / dr/e*w£j * ■ (19) 

J O77' gfrV 1 ) +? - 1 



2 Functions ^(t) and 6(f) 



Coming to the new variables in Boltzmann approximation of the function $(2;) 
(|14|) the expression for the function ^(77), in which the statistic factor is already 
taken into consideration, is 

00 

,(,)_ 4^! — (20, 

Let us introduce new variables x and z 

i = ^/rjsmh(x), z = cr^/rj. (21) 

Then we get 

7V X 3 r s i n h 2 tdt 

^' ~ J e * COBht -l' 



Calculating the integral J ^dt and changing the order of integration in the 
obtained expression we arrive at 

00 00 00 

B(r 1 )=[mdt'= 2 -^[^d X [ ^-dv. (22) 
J it (7 J cosh x J e v - 1 

t y/V " cos h ^ 
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In particular integrating over the whole interval of the value t we obtain the 
integrals product 



9(0) = / ^{t')dt 



2N X [ sinlr x 



7T cH J cosh 



dx 



-dv, 



one of them is expressed in terms of Q - Riemann function 



1 15' 



and the other is easily calculated 



sinh x 7r 
—dx = — . 
x 16 



cosh 5 



As a result we find 



9(0) = / V(t')dt 



TT 3 N X 

120a 2 ' 



Thus we can write down 



9(7?) = 9(0) - 



2N X f sinh x 



^/rjcr cosh x 



7T <T J COSh X 



dx 



-dv. 



The inner integral can be expressed in terms of the function^] 

X 

o 

the function D{x) has the following asymptotics 



D(x) 



3E 



(n + 3)n! ' 



< 1; 



F-3 -3 ( 1 + - + -a e ' x>1 ' 
oar V a; ar 



(23) 



(24) 



where £?„ are Bernoulli numbers. The function D(X) graph is shown in Fig. [TJ 



2 The function D(x) is connected with the Debye functions (Ref., for example [3])- 
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Fig.l. The function D(x) graph constructed with the help of the Maple pack- 
age on the approximations of the type $24\) . 



Thus finally we obtain 



9(77) 



s(vV)> 



120cr 2 

where a monotonically decreasing function S(x) is introduced 



E(x) = 1 



80 



■|ta n h>*D(xco B h*)* ) 



(25) 



(26) 



varying in the interval 



> 3(x) < 1. 



To calculate with the functions D{x) and S(x) a special library in the package 
of symbol Maple Mathematics was worked out in which the proceeders of time- 
optimal calculation of these functions with the help of different approximations 
are determined. The function S(x) graph is shown in Fig. [5] 
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Fig. 2. The function S(x) graph. The dotted line denotes the extrapolating 
function F(x) graph if27|) . 



In the interval [0, 10] the function S(x) is well extrapolated by the function 



F{x) 



1 + fix 2 

with the parameters a = 0, 05 and (3 = 0, 09. 



(27) 



3 Equilibrium deviation 5f(rj, £) 

Now let us calculate the function 6f(r), £). As it is not difficult to see the function 
$(77, £) is a slowly varying function, because in consequence (|15|) 



while 



K < 1/2, 



lim e* ( "' £) = 1, 



and in wide limits of changing the variables n, £ : 

cxp($(77,0) « 1. 

And the derivative of the equilibrium function of distribution has the order 



/o 2v / vTTFe <J v / ^-l 



(28) 
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and infinitely grows in the range a •*/ r\ + £ 2 — > 0; By large values of this argument 
this variable is small. Thus, values of the variables r), £ in the range \fr\ + £ 2 < 
er -1 in which exp($) can be considered approximately constant make the main 
contribution to the distribution function deviation from equilibrium. Thus, 
integrating by parts in (I19|) we hnd approximately 



Sf(r),0 



e°i - 1 



1 



(29) 



Substituting this function into the inner integral (|17[) and introducing the func- 
tion 



m(r,,a)=a 



we obtain 



u 

OO 

2vr 6 AA 



1 



- 1 



a I dr/e- e( " ) ^Df(?7,a). 



(30) 



(31) 



Note, the introduced earlier function Df(r;, a) is proportional to the perturbed 
track of tensor of X - bosons energy-impulse 



ST X = g ik 5Tg 



'x 



SfxdP. 



P(X) 



Coming to the numerical integration in (|3"Tj) let us note that the inconvenient 
for a numerical integration integrals are just in in the function Df (77, er), so the 
direct application of numerical integration runs against the divergence problem. 
Therefore, at first the integral (|3"0|) is necessary to be transformed into a conve- 
nient for a numerical integration type. For this let us rewrite the integral (|30[) 
in the equivalent form 



m(r ll <7) = a 2 




Let us study the first part of the integral 

00 

edi i 1 
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in the first integral we make a substitution £ = y/rjx, and in the second £ = 
y/x 2 — 1. Then we get 



A =77 



x dx 



VT+a? - 1 



y/ri + x 2 (£r 



Bringing these integrals together we obtain 
l 



A = 77 



x dx 



1 



Now we transform the -B part of the integral 

oo 

£ 2 d£ 1 - e-*fo'« 



£ = 



Vr? + £ 2 e^V^ 



1 



Substituting the expression for $(17, £) into the integral and making a substitu- 
tion £ = y/r]x 1 we transform this integral 



B = T) 



2 dx 1 — 



vTT^ 



In this integral the difficulties of numerical integration arise under the conditions 
of o^fr\ — > 0. Expanding the exponent into Tailor series small values 77, we get 
approximately 



1 9 
~ 2 



x 2 dx VTT^-x 2 In !+^+^ 



Numerical integration of these expression does not run to any difficulties already. 
Considering the above comments in the Maple package there was created a 
library of special procedures of rapid calculation of the function Df (77, a) for 
any values of variables. The graphs of the functions Df(?7, a) obtained with the 
help of these procedures are shown in Fig. [3j 
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Fig. 3. Functions Df (rj,a) subjecting to the parameter a. Along the abscissa 
axis lg 77 is plotted. From the left to the right a = 10, a = 3, a = 1, a = 0, 3, 

(7 = 0,1. . 



4 Results 

Before proceeding to the results presentation let us carry out a convenient func- 
tion 5s normalization. As noted in pQ, in the papers [!]-[6] the specific entropy 
estimation for one baryon was obtained (formula (5) [T]): 

_45C(3)iVx , . 

d s- 4 ^4 N Ar - i 32 ) 

Therefore we shall correspond our results to this estimation by introducing a 
relative variable 

A S - I (33, 

so called the reduced specific entropy. Carrying out numerical integration in 
the expression (|3Tj) with the help of the given procedures in the Maple package 
we obtain a graphs family of the function As (a). In Fig. 0]we show calculated 
graphs of dependence As (a) by different values of Nx - number of X— boson 
types which is a parameter of the model of interactions. 
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Fig. 4. The given specific entropy for one baryon, lg As(cr, Nx), subjecting to 
the number of X-bosons, Nx- Along abscissa axis the values of lg a are plotted. 
In the right part of the figure from bottom up Nx — 1; Nx — 3; Nx = 10; 
N x = 53. . 



Passing on the analysis of the results we firstly note that the calculations 
carried out in terms of complete kinetic theory Showed a sufficient dependence 
of the produced baryon charge in the quantity of X - bosons types. Let us point 
out the following general tendency of this dependence of the produced baryon 
charge: at er < 0,4-^0,8 with the increase of the number of X - bosons types 
the given specific entropy increases , and at a > 0,4 0,8 on the contrary it 
decreases, moreover in the range of small values of the parameter a the depen- 
dence of the given entropy on Nx is especially perceptible. At the same time 
we should remember that the absolute value of specific entropy equals to 

' S ^ s ^ 4s ffl|A, (34) 

On the other hand we can suppose the factor Nx/N (relation of the number 
of X-bosons types to the general number of particles types) does not strongly 
depend on the field model of interactions, therefore the conclusion concerning 
dependence of the number of X - bosons types of the given entropy can be 
carefully transferred to the absolute value of specific entropy also. These pecu- 
liarities of dependence of specific entropy on the number of X - bosons types 
are shown in Fig. [5] 
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Fig. 5. Dependence of the given specific entropy for one baryon on the number 
of X - bosons types. Bottom-up a = 0,1; a = 0,4; a = 0, 7; a = 1; a = 3; 
a = 10. Along the abscissa axis the values oflgNx are plotted, along the 
ordinate axis the values of lg Ag are.. 

Further, the value (|32|) obtained by a number of authors, in the kinetic 
theory is reached at the values of the parameter a = 0,4 0,6. Moreover 
the kinetic model of cosmological baryogenesis detected a finer structure of this 
process than the one which was produced by the hydrodynamic theory of this 
process developed earlier in the quoted papers [7]- [9]. The difference of our 
results from the quoted papers results, especially in the range of small values of 
the parameter a, is caused by essential influence of nonequilibrium processes on 
the final result in this range. It is not difficult to see that it is this range where 
the function of X-bosons distribution most differs from the equilibrium one. Of 
course this fact cannot be taken into consideration in the hydrodynamic model 
of bariogenesis effectively. 
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